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Abstract 

Some experimental results of correlation functions in Bose-Einstein (B-E) interfer- 
ometry measurements exhibit a non smooth behaviour - oscillations. Possible origin 
of such a behaviour in non-trivial spatial distribution of the source is considered here. 
A method for obtaining additional information about the structure of source emitting 
identical particles is suggested. It is shown that this information could be extracted 
by a one dimensional fourier analysis of the filtered correlation function. Further pos- 
sible enhancement of the method is sketched. Some technical aspects of the proposed 
technique are discussed. 

Parti: Introduction 

The idea of measuring properties of a source of identical particles by correlation experi- 
ments originated in radioastronomy. The angular diameter of the star Sirius was successfully 
determined by HBT method [1] in agreement with predictions. In experiments on antiproton- 
proton annihilation in 1959 the correlation effect in like-sign two pion angular distribution 
was discovered. This effect was interpreted in [2] as a result of B-E interference of the wave 
functions of emitted pions. In the simplest, sufficient for our purposes form the correlation 
function for a static source can be up to the normalization factor expressed as: 

C{Ap)c^l I fixi)fix2)\^uixi,X2,pi,P2)\^d'xi£x2 (1) 

where pi,p2 are momenta of emitted identical bosons and f{x) describes the geometrical 
distribution of the source. The squared absolute value of \E'i2 is: 

I P = ^ I e'P'^' e'P'"^' + e'P'^' e'^'"^' | ^ = 1 + cos( Af ■ A^ (2) 

and it depends only on the relative momentum of the emitted bosons Ap = pi —p2 and on the 
relative position of the emission points Ax = Xi—X2- For a Gaussian source f(x) ~ e~^^ 1^ ^ 
the correlation function (1) can be calculated analytically: 

C{Ap)- J J \^u\^e-^'^'-'^'^/^'£xid'^X2 = l + e-\P'-P'\'^'/^ (3) 

By fitting experimental data for the correlation function with (3) an approximate radius 
of the spherical source can be determined. A more complicated Gaussian parametrization 



of the correlation function [3] separates the longitudinal, outward and sideward dimensions 
of the source. Statistical errors of correlation function data in heavy ion collisions (HIC) 
experiments became smaller but the interesting structure of the correlation function in the 
higher relative momenta region is still present in the recent results [4]. The influence of 
Coulomb final state interaction leads to the more interesting shape of the theoretical corre- 
lation function for the like sign charged particles [5,6], but it is not rich enough to explain any 
oscillations of the experimental correlation function in the Ap > lOOMeV/c region^. Per- 
haps, it could be possible to reveal some information about the source from this behaviour 
of the correlation function. 

In the following section we define the longitudinal correlation function which is subse- 
quently used in the definition of distance spectrum of the source in Sec. 3. In Sec. 4 we discuss 
the reconstruction of the distance spectrum function from the longitudinal correlation func- 
tion and some technical aspects of this procedure. In Sec. 5. we present a simple example 
as the illustration of the method. In the last part we consider the possibility to apply a 
tomography techniques based on the approach presented here. 

In the whole paper we deal only with the influence of the spatial distribution of source 
on the correlation function. The effects of the final state interactions and the coherence of 
emission are not considered here. It is assumed that they do not affect the principles of the 
method presented. 

Part2: Longitudinal Correlation Function 

Let us consider a spatial static three-dimensional source emitting pairs of boson particles. 
Because of B-E interference effect, the emission amplitude of a pair of particles with relative 
momentum Ap = P2 — P2 depends on the relative position of the emission points. We can 
decompose the relative momentum of the detected pair into transverse and longitudinal 
components Ap = Apt + Api, where the longitudinal direction is chosen along the beam in 
the HIC experiment. Thus we have: 

Ap-Ax= (Api + Apt) ■ {Axi + Aft) = Api ■ Axi + Apt ■ Axt (4) 

Let us choose events with a small transverse component of relative momentum 

\Apt\ « \Api\ (5) 

Then the transverse separation of emission points does not have influence on the corre- 
lation function for our set of data. An analogous consideration can be made for the set of 
events filtered by the condition 

\Apt\ » \Api\ (6) 

or one can choose the sideward and outward selections used e.g. in [3] but we shall now 
concentrate on the correlation function for data events fulfilling the condition (5) - the 
longitudinal correlation function. In this case after the substitution (4) into (1) and assuming 
(5) we have: 

C{Ap) = J J f{xi)f {x2){l + cos{Ax- Ap))d^xid^X2 = 

= J J F{xj)F{xf){l + cos{Axi- Api))dxjdx^ (7) 

^ The oscillations in theoretical shapes of the correlation functions for photons can be found in the paper 
[7] (but see also [9]) 
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where F(xi) = J f(xi,Xt)(fxt. 

This means, that our 3D source produces C(Ap;) hke a one dimensional source with the 
distribution F{xi). 

PartS: Distance Spectrum of the Source 

Let us find how the longitudinal correlation function C{Api) looks like for a two-point 
source (see Fig. la). This two point source has a distribution function F{x) = l/2{5{x — a) + 
6{x — b)). The longitudinal correlation function then is: 

C2{d, Api) = j F{xi)F{x2){l + cos{Api{xi - X2)))dxidx2 = 1 + cos(Ap, • d) (8) 

where d = \a — b\ is the separation of emission points. For a source with the distribution 
function F(x) = Pi ■ 6{x — xi) + P2 ■ 6{x — X2) + P3 ■ 6{x — xs) (Fig. lb), where Pi + P2 + P3 = 1 
the result is: 

C{Api) = e + P12 • C2{du: Api) + Pi3 • C2((ii3, Api) + P23 • C2((i23, Api) (9) 

where Pij — Pi-Pj, dij — \xi — Xj\ and © is a constant shift which will not be important in our 
further calculations. Based on this it seems appropriate to define the Distance Spectrum of a 
source. Each source emitting pairs of particles has its 'squeeze' - Distance Spectrum. It is the 
probability distribution for the emission of the pair in the distance D of emission points. For 
the two-point source the distance spectrum is S{D) ~ (5(P> — do) (see Fig. la), for our three- 
point source (8) the distance spectrum is S{D) ~ Pi2S{D—di2)+Pi35{D—di3)+P235{D—d23) 
(see Fig. lb). In general the distance spectrum for a one dimensional static source can be 
expressed as: 

SiD) = I f{xi)f{x2)5{D - \X-, - X2\)dXidX2 (10) 

Now we can write the longitudinal correlation function produced by a source with distance 
spectrum S{D): 

C{Api) = Q' + J S{D)C2{D, Api)dD (11) 

A distance spectrum does not contain complete information about the spatial distribution 
of a source. Slightly different spatially distributed sources can have the same S{D) (see also 
Sec. 5.). The distance spectrum of another point-like source is shown in Fig.lc . 

PcU"t4: Inverse Transformation 

Expression (1) for the correlation function can be after the formal integration written in 
the form: 

C(Ap-) = l + |/(Ap-)|2 (12) 

where f{Ap) = J e^^^'^ f{x)d^x. The absolute value in expression (12) destroys the phase 
information in the fourier picture f{Ap) of source distribution. This breaks the possibility 
to perform the inverse fourier transformation in order to get f{x). 

In this section we shall show that the remaining information in |/(Ap)| has certain 
physical meaning and that its extraction is at least theoretically possible. 
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The form of expression (11) is sufficiently simple. Therefore one could consider a possibil- 
ity to gain the function S{D) from C{Api). As we shall see, the S{D) is just the information 
contained in C{Api). According to (11) and (8) we have 

C{Api) = e" + J S{D) cos{Api ■ D)dD (13) 

The constant shift factor can be separated from the correlation function data and so we 
can write: _ 

C{Api) = / S{D) cos{Api ■ D)dD (14) 

JO 

where C{Api) is the correlation function with the shift factor removed. Expression (14) is 
fourier transformation and we can try to perform the inverse fourier transformation in order 
to obtain S{D) from experimental data. 

S{D) ~ / C{Api) cos{Api ■ D)d{Api) (15) 
Jo 

Because of the fact, that we measure C{Api) only in discrete points it is necessary to 
replace integration by summation. 

points 

SpiD) ^ Y: C{Api) cos(Apj ■ D){Api-' - Ap]) (16) 

i 

However a more relevant problem is that the integration region in (15) includes the 
values of Api which cannot be measured ( high Api ). The cut-off in Api is present due to 
energy conservation considerations and due to the growth of statistical errors in a high Api 
region. The influence of such a cut-off on the result of Fourier transformation leads to the 
oscillations in the resulting S{D) (see Fig. 2a) regardless of the number of measured points in 
the region {0,Ap^^^). Similar problems occur e.g. in digital FIRQ filter design. Fortunately 
these oscillations can be suppressed by the Method of Windows [8] . The method is based on 
the simple multiplication of the original function to be fourier transformed by the Window 
Function, which suppresses the values of the original function near the cut-off (see Fig. 2b). 

points 

Sw{D) ^ Y WiApi)CiApi) cos(Apj ■ D)(Apr' - Arf) (17) 

i 

There are several types of the Window Function. In Fig. 2b a simple Gaussian function 
is used. The amplitude of the oscillations depends also on the shape of the exact (without 
cut-off) result of the Fourier transformation. The exact result of Fourier transformation of 
the function in Fig. 2a without cut-off would be Dirac delta function at the point rfo = 3fm. 
Such a shape leads to a big oscillations. 

The influence of statistical errors of correlation function data to the resulting distance 
spectrum which is crucial for the eligibility of the method is not considered in this paper. 

Part5: Structured Time Dependent Source 

As we have already mentioned the sources with the same S{D) can be different, but in 
spite of that, having the distance spectrum of a source we can make some conclusions about 

^FIR — finite duration impulse response 
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the source. At least the average radius 

R = J D- S{D)dD (18) 

and maximal radius i?max = -Dmax in S{D) can be found from S(D). 

Let us consider a three-dimensional source consisting of two separated regions (Fig. 3c). 
Both regions A and B contribute to distance spectrum of a source in longitudinal (5) analysis 
in the interval < D < R. When one of the particles comes from A and the second 
from B they can be emitted in a distance interval K — R<D<K + R. The distance 
spectrum S{D) in the interval R < D < K — R is suppressed as much as the parts A 
and B are located. Original spectra Smc{D) for such a source computed by Monte-Carlo 
simulation, the correlation functions calculated according expression (11) at 200 points in the 
Api < 2GeV interval and the results of inverse transformation (16) for different separations 
K of parts A, B of the source are in figures 3a, 3b, 3c. 

The real experimental situation is however a time dependent distribution of the source. 
The general expression for the correlation function in this case is much more complicated 
than expression (1): 

C{Ap) = J J w{x^, Pl±Pl)w{x,, ^^^)(1 + cos{Ax^Ap^))d^x,d'x2 (19) 

However the influence of the spatial structure of the source could still demonstrate itself 
in the oscillations of the correlation function at the higher relative momentum region. 

Following the simple approach based on (1) in the case of the time evolution we have to 
take into consideration that the bosons emitted in different times do interfere. This means, 
that the same region of the hadron gas or quark gluon plasma emitting the bosons in two 
different times ti and t2 at the positions Xi and X2 contributes to our distance spectrum of 
the source hke a two different regions emitting the bosons at the same time but in separate 
points xi and X2- Therefore an important question arises: How big can be the time difference 
of the processes of emission for interfering bosons ? 

Let us denote the duration of the process of emission of a single boson as Tq. As an 
upper limit for tq the formation time of the emitted boson which is of the order of Ifm/c [11] 
could be chosen. It seems to be reasonable to require a time overlapping of the processes 
of emission for the interfering bosons in the B-E interference phenomenon. Assuming the 
Lorcntz dilatation factor 7 for tq of emitted boson we obtain the rough restriction on the 
time difference At of emissions of the interfering bosons in the form: 

At < To • 7 (20) 

For a typical momentum of pions produced in heavy ion collisions the time At is long 
compared to what is necessary for the influence of B-E interference on correlation 
function [10]. 

Different spatial distributions in a different time intervals lead to different distance spectra 
s{D,t). In our simple approach the correlation function for a time dependent source could 
be expressed for tq — > as: 



_ r rco poo [ n 

C{Api)^JJ s{D,t)cos{Api- D)dDdt^ j y s{D,t)dt 



cos{Api ■ D)dD (21) 



but for To 7^ the information about the s{D, t) in the correlation function is smudged with 
the uncertainty proportional to At given by (20). 
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It is clear, that by the inverse transformation (15) we can obtain only S{D), the result 
of the time integration in (21). A combination of the results obtained for different kinds of 
identical particles emitted in different stages of the evolution of the heavy ion collision or 
the results obtained for filtered events which can be produced only during characteristic and 
short time intervals of the evolution process could bring some information about s{D,t). 

Part6: B-E Tomography 

In this section we shall deal with the set of events fulfilling condition (6) - transversal 
events. Analogously to Sec. 2. the longitudinal separation of emission points does not have 
influence on the transversal correlation function^ now, but condition (6) is not so selective 
as condition (5). We have a freedom in azimuthal direction of Apt and the emission points 
lie in the plane orthogonal to the beam direction (see Fig.4). It is possible to apply some 
further condition to have the same one-dimensional situation like in the preceding sections. 
For example the condition 

Apt ■ » \Apt X (22) 

where normalized vector orthogonal to the beam direction (see Fig.4), selects the 

pairs with relative momentum in the direction of vector n^. 

However the cylindrical symmetry of the radiating volume in HIC, at least for the central 
collisions, instigates us to try to gain the whole two-dimensional information from C{Apt). 
Let us define two-dimensional distance spectrum as: 

S{D) = J J f{xi)f{x2)6{D - (xi - X2))d^xid^X2 (23) 

S{D) is the probability distribution of emission of the pair of particles from the points at 
relative position xi — X2 = D on the transversal plane. Then for our transversal correlation 
function 

C{Apt) = J f{xi)f{x2) cos(Ap4 ■ Ax)d^xid^X2 (24) 

we can write 

C{Apt) = J S{D) ■ cos{Apt ■ D)d^D (25) 

what can be proven by inserting (23) into (25). 

For a set of transversal events fulfilling the condition (22) specified by the angle (p the 
correlation function is: 



C{Ap^) = / / S{D) cos(A^^ ■ {D^ + D^)dD^dD- 



S{D)dD, 



_L 



cos(Ap^ ■ D^)dD^ = (26) 



J S^{D^) ■ cos(Ap;p ■ D^)dD^ 



where the index " describes the orthogonal and the index " ^" parallel direction to the 
vector n. Thus by the inverse transformation (15) we can obtain the projection S^p{D) of the 
two-dimensional distance spectrum S{D) to the direction determined by the angle ip. This 
can be done for any direction, for example in 2° steps of our angle (p. 



the correlation function for the transversal set of events 
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From such a set of projections one can get the whole S{D) by Radon Transformation 
which is used in tomography [13]. 

It is hard to imagine that this technique could be applied to the interferometry exper- 
imental data in heavy ion collisions. However in the future HlC experiments e.g. at LHC 
the multiciplities/event could allow to do something close to the ideas presented here. For 
the final decision whether this tomography method is eligible or not a careful and precise 
estimate of statistical errors is necessary. Both the subjects - the description of Radon 
Transformation and the analysis of statistical errors are not included in this paper. 

One can imagine that the proposed two-dimensional tomographic method is applicable 
and even more, that a time dependent three-dimensional distance spectrum can be somehow 
gained. The question arises: What could we see from S{D) or S{D,t) ? Author thinks that 
it could be possible to see some signatures of the phase transition in this kind of analysis. 
In cosmology the phase transition as the only explanation of the inhomogeneities in the 
distribution of matter in universe is used. We can hope that the phase transition in HIC 
could lead to the formation of the inhomogeneities - the bubbles in the volume of the collision. 
These bubbles emitting particles in a different strength than the surrounding medium could 
be seen in our distance spectra. 

Part7: Summary 

We have studied the influence of the spatial distribution of a source on the correlation 
function in a set of longitudinal (5) events. It was found, that spatially complicated structure 
of the source can lead to oscillations of the correlation function. It is shown that from the lon- 
gitudinal correlation function the distance spectrum of a source in the longitudinal direction 
can be derived. The radius of the source, which is used to characterize the longitudinal size 
of the collision volume is expressed as a simple integral of the distance spectrum. Some tech- 
nical aspects of the suggested method are discussed. The problem of the time dependence of 
spatial distribution of a source is also considered. An analogous method for transversal set 
of events and the possibility to use the tomography method in this approach are considered. 
The proposed technique could evoke some ideas how to enhance the interferometry methods 
in the future heavy ion collisions experiments. 
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Figure captions 



[Fig. la,lb,lc] The point-like sources and their distance spectra obtained by Monte-Carlo 
simulation in the agreement with the analytical probability calculation. 

[Fig. 2a, 2b] The correlation function C2{Ap) for two-point source (distance = 3fm) and 
the result of the inverse transformation without (2a) and with (2b) the Method of 
windows. 

[Fig. 3a,3b,3c] The distance spectrum Smc{D) calculated by Monte-Carlo simulation, the 
corresponding correlation function C2(Ap) at 200 points and a result of inverse trans- 
formation (15) for the three different separations of the sources A, B. 

[Fig. 4] The plane transversal to the beam direction and the vector n with direction de- 
termined by the angle (p. For the central collision the zero angle direction can be 
chosen randomly. For a non central collision an asymmetry in azimuthal direction in 
the distributions of the produced particles [12] could be used to determine the zero 
angle direction. 
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